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NOTATIONS 


a length, of rectangular plate 

b width, of rectangular plate 

h thickness of rectangular plate 

x coordinate along plate length 

y coordinate along plate width 

a x ,cTy normal stress components 0 X and a y generally counted posi- 

tive when compressive stresses 

t shear stress component 

e x ,£y normal strain components e x and e y generally counted posi- 

tive for compressive strains 

7 shearing strain component 

a critical state of stress 

6 critical state of deformation 

a e critical stressed state computed by elasticity theory 

e e critical state of deformation computed by elasticity theory 

3 X critical stress 0 X of a rectangular plate under uniform load 

in the length direction 

i x corresponding strain e x 


"Over het knikvraagstuk in het plastische gebied bij staven en platen 
(Deel II)." Nationaal Luchtvaartlahoratorium, Amsterdam, Rapport S. 423, 
Jan. 1954. 
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5 X e ,€ x>e the same, computed .by elasticity theory 

critical shear stress t of rectangular plate loaded in shear 
corresponding shearing strain y 
the same, computed by elasticity theory 


T 

7 

*e>h 

k 


so-called buckling factor — a nondimensional factor appearing 
in one formula in the form (see section 2.2) 


a = 


kit 2 Eb 2 


12(l - v 2 )h 

buckling factor computed by elasticity theory 


t) = a / 5 e = k / k e 


E 


E r 


Young's modulus — 

tangent modulus , that is, dff x /de x of the stress-strain 
diagram of the tension and compression test 

secant modulus, that is, cr x /e x of the stress-strain diagram 
from the tension and compression test 

reduced modulus according to (3-5) 


a E/E r 

G -shear modulus 

v Poisson's ratio, that is, -e x /£y of the bar stressed in 

tension or compression along x into the elastic range 

v' the quotient of -€y/e x of the bar stressed in tension or 

compression along x into the plastic range 


2. INTRODUCTION 


Part I of this report (ref. 1) dealt with the buckling of plates 
in the elastic range on the basis of two different kinds of plasticity 
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theories, the so-called deformation theory and the so-called incremental 
theory. 

It is now generally accepted that, from the physical point of view, 
incremental theories give a rather good picture of the plastic deforma- 
tion (compare refs. 2, 3, and 4) and that deformation theories certainly 
cannot be correct (refs. 2 and 3)* However, it still seems that the 
results for the critical buckling loads of plates, loaded in the plastic 
range, are substantially higher when computed by incremental theory than 
indicated by the test . 

The experimental values are correct in the neighborhood of the 
results by deformation theory, although quite often the "excellent 
agreement" spoken of does not exist. 

The failure of incremental theories for defining the critical buckling 
load in the plastic range is frequently regarded as proof of the inaccuracy 
of incremental theories themselves, but that is entirely unjust . As Prager 
so rightly remarked (in ref. 5): "To this writer (Prager) the idea of 

testing a stress-strain law by buckling experiments seems utterly fan- 
tastic; nobody would dream of determining, say. Young's modulus by a 
buckling test in the elastic ' range rather than a simple tension test. 

Direct experiments speak in favor of incremental theories, and there are 
strong theoretical objections against deformation theories." 

Batdorf and Budiansky (ref. 6, see also refs. 3 and 1) developed a 
theory which, for certain load increases, gives results according to a 
deformation theory. Their theory has less physical faults than deforma- 
tion theory itself, but some direct checks on this theory come out badly. 

Besseling (ref. 1) sums up some of the possible causes of the fail- 
ure of incremental theories and also cites the report by Onat and Drucker 
(ref. 7). The last-named report shows, on a special load case, how small 
Initial eccentricities exert a very strong reducing effect on a critical 
buckling load computed by incremental theories, as shall be explained in 
section 3«5* That section deals also with a special investigation that 
indicates how the stress-strain relations by both theories vary when, 
at a constant compressive load above the elastic limit, a shear stress, 
increasing from zero, is applied. 

In practice, it is hardly possible, of course, to include Initial 
eccentricities in the calculations; they should also be exactly known. 

The fact that then all points of a plate have different stress- strain 
relations (which in the special load case used by Onat and Drucker was 
not the case) raises well-nigh prohibitive mathematical difficulties. 

The question of whether the calculation by deformation theory is the 
right way forms the subject of the present report. Yet, before applying 
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deformation theories to buckling problems of plates — and this is accom- 
plished before passing to incremental theories — some simple calculating 
methods are given which, compared with the experiment, did not give such 
bad results. 


2. STATEMENT OP RESUITS 


2.1 INTRODUCTION 


On a rectangular plate, loaded in its plane by a uniform stressed 
state, the components of which are defined by a constant factor, the 
critical state of stress indicated by 5 is defined by length a, 
width b, plate thickness h, by the clamping at the edges and by the 
material properties. Variation of the linear plate dimensions in the 
same ratio leaves the stress or strain components in the critical state 
unchanged . 

By "load case" is meant a case in which a rectangular plate of 
specific material and with specified attachments at the edges, changes 
to a uniform state of stress before reaching the critical state, and 
whereby the stress components have a certain fixed relationship. 1 To 
illustrate: a rectangular plate of a certain material, clamped at the 

edges in uniform state of shear stress, is a load case. Another one is 
the rectangular plate of a certain material with hinge-supported edges 
in uniform state of normal stress (with respect to a pair of sides) . 


2.2 Variation of Length-Width Ratio 

The critical stressed state 5 for a specific load case is, as far 
as the plate dimensions are concerned, simply dependent on two independent 
variables, such as the length-width ratio a/b or the width-thickness 
ratio b/h. 



Figure 2.1 represents a diagram for constant width-thickness ratio, the 
critical load plotted against the length-width ratio. 


Only such load cases are treated. Plates other than rectangular 
are excluded, so also plates whose stressed state before reaching the 
critical stressed state is not uniform, although defined by a constant 
factor, such as the web of a beam under moment load or shear, for example. 
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The different branches of the curve correspond to different num- 
bers (m) of half-waves. The number of half-waves computed by different 
plasticity theories or elasticity theory need not be identical for a 
given load case a/b. 

By (3*1) it follows that the critical stress 5 for a given load 
case in the elastic range can be written 


5 = jj- g( a jb) = ^ j(a,b) 


( 2 . 2 ) 


where g and $ are functions of the rectangular sides a and b. 
Examination of the dimensions of a, D, and h indicates that func- 
tion j is dimensionless. Thence by (3-2) 


- _ rt 2 kP _ kjt 2 Eh 2 
kh 2 12(l - v 2 )b 2 


(2.3) 


with k a function of the length-width ratio a/b. The formula (2.3) 
is, naturally, in agreement with the form of (2.1). 


Formula 


3 = 


ha 2 


with 


k' 



(2.4) 


is, of course, valid too. 

The square of the side b loaded in compression is chosen in the 
denominator of formula (2.3) for the plate compressed in one direction, 
solely because otherwise the limit transition to a— becomes 
difficult . 

The coefficient k in (2.3) is called the buckling factor. 

The critical stress state in the plastic range is indicated by a 
different function of the ratio b/h than given in (2.3) , though the 
form of (2.3) is also used in this range, but then the buckling factor k 
* is no longer a function of a/b but of b/h too. 
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2.3. VARIATION OF THE WIDTH-THICKNESS RATIO 
2.3 .1. Kbllbrunner Method 


When setting up the length-width ratio, the width-thickness ratio 
remains merely as variable. The length-width ratio of an infinitely long 
plate can equally be considered as fixed when the critical load at 
approaching » plate length does not approach zero. 

Figure 2.2, taken from Kollbrunner 1 s report (ref. 11), represents 
the critical load of an infinitely long plate supported by hinges on the 
long edges and under lengthwise compression plotted against what 
Kollbrunner calls the slenderness 


A = £0-2(1 - v 2 ) 


2.3.2. U. S. Methods 

2.3.2.I .- When the length-width ratio of a plate is fixed (for the 
plate of infinite length this means, if the condition cited in section 2.3 
is complied with), the NACA does not select the width -thickness ratio b/h 
as independent variable, but the critical state of stress or critical state 
of strain, as computed by elasticity theory. The thus computed states are, 
for a certain load case, with fixed length -width ratio, naturally simple 
functions of the b/h ratio, and therefore can supersede the last named 
ratio as variable. 

In figure 2.3 the (actual) critical state of stress (vertical) is 
shown plotted against the critical state of strain computed by elasticity 
theory (horizontal). When, as in figure 2.3, the critical strain condi- 
tion computed by elasticity theory is plotted horizontal, the NACA us ually 
employs the term "strain" or "shear." This simple indication is often too 
Indefinite, since on the vertical the term "stress" or "shear stress" is 
used also (i.e., the actual critical stress or shear stress). Some writers, 
like Schuette and MacDonald (ref. 10), even call the curve in a diagram 
with such indications the "plate compressive curve," by which the impres- 
sion is gained as if the diagram gave the relationship between stress and 
strain (or crushing of plate divided by height = mean strain) of- a plate 
during a loading process. - ; 

By chance the curve in a diagram of the discussed calculating method 
by Gerard (section 3-2) was exactly identical with the stress-strain curve 
of the material. The true significance of the diagram is, however, an 
entirely different one; it gives the critical state of stress in terms of 
the width -thickness ratio. Some writers who use other computing me thods 
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than Gerard's frequently plot, besides their computed curve, or besides 
test points, the actual stress-strain curve of the material. This has 
the advantage of making the value of the critical strain available while 
reading off a critical stress value. 

2.3«2.2 .- From a diagram as discussed in section 2.3«2.1 where a 
critical stress a is plotted against a e computed by elasticity theory, 
another one can be constructed in which 

r] = — against a 
5 e 

This form of diagram is also much used by the NACA. cr (horizontal) 
against d/a e (vertical) . 

In this manner figure 2.4 was derived from figure 2.3. Point A in 
figure 2.3 has the coordinates: 0.008 and 60,000 lb/inch2. Young. 1 s 

modulus is E = 10. 7 X 10 ^ lb/inch^. The corresponding point A in fig- 
ure 2.4 has then 60,000 lb/inch^ as abscissa and 

T) = 60,000/0.008 X 10.7 X 1C)6 = 0.7 as ordinate-. The factor t] can be 
expressed by 



with k the buckling factor and k g the buckling factor by elasticity 
theory. 

The factor t] is always the factor with which the critical load 
computed by elasticity theory must be multiplied in order to obtain the 
actual critical load. Formulas applicable in the elastic range cannot 
only be generalized by giving the buckling factor k another value than 
by elasticity theory, but alsq by retaining the value of k by elastic- 
ity theory when adding the factor q to the particular for mula . The 
factor tjE appearing then in such a formula is called the buckling 
modulus, as by Kollbrunner in reference 8.2 

2.3 .2.3 Besseling (in ref. 1) uses a diagram in still another 
form; not t] = d/a e plotted against 5, but buckling factor k = Tjkg 
plotted against critical strain I, which, however, is associated with 
the critical stress a in very simple manner. For the plate confessed 
in one direction it naturally gives e = 5/E sec . Figure 3*4 represents 
such a diagram. _ 

^he same writer calls the tangent modulus "variable elasticity 
modulus" in reference 8. 
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3. THEORIES OF PLASTICITY USED TO DEFINE THE CRITICAL 
BUCKLING LOAD OF PLATES 


3.1. INTRODUCTION - SHANLEY EFFECT 


If a plate, loaded in its plane only, and hence "being initially in 
the state of stress of a membrane, is subjected to a deflection perpendi- 
cular to its plane, bending and torsion stresses arise, and therefore 
the stress or the criterion for further plastic deformation (generally 
the distortion energy) over a section of the plate thickness decreases. 

If the plate in its uniform state of stress was already loaded as far as 
in the plastic range, the appearance of a deflection perpendicular to the 
plane of the plate makes the corresponding deformation partly plastic and 
partly elastic, which complicates the calculations. 

Shanley (ref. 12), however, pointed out that the first appearance of 
a deflection in consequence of reaching the critical buckling load is 
still accompanied by further simultaneous deformation in the plane of 
the plate, so that progressive plastic deformation takes place everywhere. 
In fact, there is no buckling process at all in the ordinary sense of the 
word (i.e., a process that takes place while the outside loads remain 
constant ) . 

Shanley applied this principle to the straight bar under column load, 
and established that the start of the buckling process can be identified 
by replacing Young's modulus E in the formulas for the buckling load in 
the elastic range by the tangent modulus E^ . As a rule, any further 
deflection is accompanied by a slight load increase only, so that the 
load at which deflection can start may be defined as "buckling load." 

This had been done In the past, obviously unconsciously, but later on 
E r (a reduced modulus) lying between E and Ej^ was used, as the 
buckling process was considered a constant load process in which, in 
fact, unloading takes place in a part of the bar. __ 

. This calculation with "Shanley effect" is found also In the treat- 
ment of plates. Nowadays It Is assumed that the calculation can be car- 
ried out as if the Shanley effect over the whole plate did occur. In 
other words, a critical buckling load In the plastic range is calculated 
solely by stress-strain relations applicable to progressive plastic defor- 
mation. But the result will not be such, even by approximation, that E 
can be replaced by E^ In the elastic formulas, in order to get the 
desired result as Is the case in the buckling of bars. 
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3.2. FIRST ATTEMPTS 


The differential equation of a buckled isotropic plate loaded in 
its plane has, in the elastic range, the form 


JzK 0 

D _ + 2 

W 4 


d 4 w 

dx 2 dy 2 


d 2 w 


bK\ . 0 ^ 

+ — r + ffy h — + 2-rh 

dy 4 / dx 2 dx dy 


+ oyh 


d^ 

dy 2 


= 0 


(3.1) 


where 


D = 


Eh-' 


12(1 - v 2 ) 


(3-2) 


E Young's modulus 

h plate thickness 

v Poisson's constant 


w displacement perpendicular to plate surface 

o x ,0 normal stresses along rectangular sides (positive, when com- 

pressive stresses) 


t shear stress along rectangular sides (just like cr x and dy 

the positive direction for t is opposite to that of nor- 
mal convention) 

Bleich (ref. 13) gives for the plate compressed in one direction 
( 0 y=O, t=0), loaded in the plastic range, after a kind of intuitive 

reasoning 


in which 



d w 

i? 


+ 


d 4 w d 4 w\ 
dx 2 dj^ Sy 4 / 




° xh d>T 2 


= 0 


(3.3) 


a = E r /E 


(3* 4 ) 


the reduced modulus 
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E r = 


4EK fc 


(Vb + 


(3.5) 


just as it applies also to the buckling of rectangular bars (when 
assuming, as done in the past, that Truckling is not coupled with pro- 
gressive plastic deformation over the entire cross section) and D 
according to formula (3*2) hence with the elastic values for E and v.^ 


Chwalla (ref. 15) assumes that E can be replaced by E r = aE 

and v by its plastic value v 1 in (3*1) and- (3.2); in other words, 
that the plate remains isotropic. He reports that the elastic value 
can also be used, for v < v' < l/2, and the 'variation in 1 - v 2 does 
not amount to much. 


When v = v 1 it naturally follows that 


= a 0. 


x,e 


(3-6) 


where a x is the critical value of cr x for the appearance of buckling, 
and 0 the value as computed by elasticity theory. In consequence. 
Young's modulus E can be replaced by the reduced modulus E r according 
to (3-5) in the formulas for the elastic range. 


Timoshenko discussed the points of view of Bleich and Chwalla in 
reference 14 . 


The concepts of Bleich and Chwalla were not based upon specific laws 
of plastic deformation, but were obviously intended as a rough intuitive 
and safe approximation. Still, the reprinting of references 13, 14, 
and 15 following the publication of Nadai's book (ref. l6), with its 
simplest form of deformation theory for a hardened material, which, 
with all Its defects, nevertheless proceeds from explicit assumptions, 
in conjunction with other findings (section 3-3) resulted in a more 
rigorously developed theory for the calculation of the critical buckling 
load of plates. 


^Many writers take Poisson's transverse contraction coefficient v, 
so that it varies in the plastic range, and at progressive plastic defor- 
mation approaches l/2 (since plastic volume change cannot occur) . How- 
ever, it is desirable to restrict v to the elastic part of the defor- 
mation, so that both v and E remain constant. The first -mentioned 
transverse contraction coefficient which thus refers to total deformation 
is expressed here by v'. 
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In Stanley's latest point of view, Bleich, Chwalla, and Timoshenko 
should not have employed the reduced modulus E r hut the tangent 
modulus E^. ' Thus, according to Chwalla, Et should replace E for the 
critical load a x in the formulas by elasticity theory. This suggestion 
was already made by Langhaar (ref. 17)* The results obtained would, in 
consequence, be too low (except -when the elastic limit itself is not yet 
exceeded), hence would give a much too unfavorable presentation of the 
subject . 

Kollbrunner ' s method (ref. 8) lies about midway between that of 
Bleich and Chwalla, as seen from his differential equation 


d^w / .1 — \ 


8^w + d^w 

dx 2 dy 2 8y* 


8x 2 


(3-7) 


Lastly,- in order to check the agreement of his tests (section 4.3) 
with the thus developed theory, he introduced a rather arbitrary correc- 
tion factor in the final critical load formulas. 


Lundquist (ref. 18) also used a formula of Chwalla (3.6) but with 
a replaced by a correction factor, while 


a + 3^ja 


This is intended as a kind of average of the aforementioned theories 

Gerard (ref. 19) substituted the secant modulus E s for E in the 
formulas by elasticity theory, without, however, giving a more or less 
theoretical basis like the aforementioned writers. According to Gerard's 
method, it can be shown that the critical state of deformation of the 
plate in a certain load case is independent of the stress-strain rela- 
tions. In figure 2.3 the curve by Gerard's method is shown as the stress 
strain curve of the material. 


3.3- DEFORMATION THEORY 


Ilyushin (ref. 20) treated the stability of plates by deformation 
theory in its simplest form (briefly indicated by deformation theory, 
because more general deformation theories are not involved here) i.e., 
with Nadai's formulas (ref. 16) where v = l/2, and with Shanley effect 
disregarded . 
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Kollbrunner, after his attempts discussed earlier (ref. 8), became 
immediately very enthusiastic and in reference 11 reexamines his own 
tests of reference 8 in the light of Ilyushin's report. He writes "The 
work of Ilyushin is theoretically correct and generally valid," but it 
was a premature conclusion. A disproval of incremental theories is not 
given j they are not even mentioned. 

Stowell (ref. 21) deals with the work of Ilyushin, still uses 
v = l/2, but includes the Shanley effect. 

Bi jlaard (ref. 22) gives the treatment of the stability of plates 
with the correct (elastic) value for v and allows for the Shanley 
effect. Although his theory has not been generally known before this 
publication, it already had been published in different form in 1938 
(in ref. 23), hence before the Shanley report (ref. 12). 

Subsequently Stowell (together with Pride) supplemented his work 
(ref. 21) with other values of v (ref. 24). With this he then got 
the same results as Bljlaard in the earlier report (refv 22) . Stowell 
compared the results of this report (ref. 2b) with those of his report 
(ref. 21); the differences seem to be very small. Bijlaard, however, 
indicates (in ref. 25) that the correction by reference 2b or reference 21 
for one of the two load cases investigated by Stowell has the opposite 
sign. 


3-4. INCREMENTAL THEORIES 


Handelman and Prager (ref. 9) applied Prager's incremental theory 
in its simplest form (in ref. 26, this theory Is called incremental 
theory, for short; more general Incremental theories are disregarded) 
to the buckling of plates without due regards to the Shanley effect. 

The exact value for v is used. 

Hopkins (ref. 27), Pearson (ref. 28, for v = l/2) and Besseling 
(ref. 1) do this, while taking the Shanley effect into . consideration. 
These last three reports should be regarded as being founded on the best 
theoretical basis. 
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3.5. DIFFERENCE BETWEEN DEFORMATION AND INCREMENTAL TEECBY 
EFFECT OF INITIAL ECCENTRICITY 
3. 5*1- Stress-Strain Relations 


After, say, a hollow tube is stressed in tension or compression in 
the plastic range with a stress Oq, an addition of torsion produces a 

distinct difference between deformation theory and incremental theory. 


If, for instance, after simple tension or compression, while the 
stress a remains constant, an increasing shear stress t is added, 
the differential quotient dj/dT by deformation theory (for proof see 
supplement (A. 3) is 


dr = 2v - 1 + 3 _ + 3g o (1 _ 3t 2 Y ,l _ _i_^ 

dT E ^ o a + 3 t 2\ a 0 2 + 5t 7 V F sec( g o) 


(3-8) 


■E's ec ( 0 o)'th‘ e secant modulus for stress cTq. 

By Incremental theory (for proof see ref. 7 or supplement (A. 2) 


dr = 2(1 + v) 9 t 2 ( 1 _ l\ 

dT E + 3^.2^ E J 


(3.9) 


Originally, i.e. at t = 0, the aforementioned differential quotients 
have the values 


Deformation theory: 


Incremental theory: 

By (3.10) the behavior is inelastic, by (3-11) elastic, hence a con- 
siderable difference between the two theories. The elastic behavior 
according to (3-H), that is to say, according to (3-9) at t = 0, does 
not continue very long, as a rule, at increasing t, for the fraction I/% 

Is generally great, and even Infinitely great for the non-hardened or no 
longer hardened material. 


dr 

dT 


2v - 1 
E 


E, 


sec 


W) 


( 3 . 10 ) 


dr _ 2(1 + v) 
dT E 


(3.U) 
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But the differential quotient dy/dT for compressed plates is 
usually related in a large measure to the critical load condition. 

Critical loads computed by incremental theory are always consider- 
ably higher than the experiment indicates. Calculations by deformation 
theory are usually in good agreement with the experiment . The failure 
of incremental theory is chiefly attributed to initial eccentricity. 
Therefore, shear stresses occur on a compressed plate before a critical 
load condition is reached. 

An explanation for the fact that experimental values for critical 
loads are in- good agreement with deformation theory (initial eccentricity 
discounted) might be found by attempting to point out that the differential 
quotient dy/dT in (3«9).> by incremental theory, already assumes, at com- 
paratively small values of t, a value equal to the value obtained by 
deformation theory at r = 0 according to (3. 10). 

If, viewed from t}ie physical point of view, deformation theory is to 
give the correct stress-strain relations, then . dy/dT should show con- 
siderable independence of t according to (3.8). Deformation theory 
gives good results even without allowing for initial eccentricity. 

Tables 3.1 to 3-8 show values for the differential quotients 


E 



and 



by incremental and deformation theory, i.e. (3*9) and (3-8). The eight 
tables apply to eight different values of Et, while v = 0.3* Each table 
shows seven values for the shearing stress r. In the choice of the values 
for ^sec^o)* °* > E(dy/d t)^ alone is dependent 


E. 


sec 


(“0) 


> E, 


must, of course, be valid. 

Definite conclusions are difficult to draw. At increasing t, 
E(d7/dT)p remains, in most cases, considerably below the corresponding 

value E(d7/d-r) D , at t = 0 (compare the numbers of the first column in 

each table with the number on top of the other columns) . For the small 
values of E^, paired with high E s g c (ff0) only, and not too low t, the 

opposite is the case (cf . table 3.? for E^ = 0.02E and compare E(d7/dT)p> 
at t = O.O 30 Q with E(d7/dr) D at E sec = 0.S® and t •= 0) . 
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For the rest it is known that 



at E 8ec (o 0 ) = E (hence at Cq = a Q , the proportional limit ) and t = 0, 

which is seen from the tables, and up to what result the other numerical 
values of the tables must converge. 


The quotient E(d 7 /dT)^ does not always increase much slower with 
increasing t than E(d 7 /d-r)-pj for high values of E sec ^ 0 Q^ the process 

is even faster, as seen, when comparing the first column in table 3-7 
with the other columns. 


3.5-2. The Buckling Case of Onat and Drucker 

From the foregoing ( 3 . 5-1) it is readily apparent that the discus- 
sion of stress-strain relations alone affords but a partial answer for 
predicting the experimental critical load by incremental theory. 

Now, Onat and Drucker (ref. 7) made some calculations on an idealized 
case of torsional buckling of a cylindrical structure consisting of flat 
plates, so that the perpendicular section represents the form of a hollow 
cruciform. (See fig. 3«1 for section, figs. 3-2 and 3-5 for bar in 
straight and twisted state.) The elastic limit for compressive stress 
is a*. Above it the tangent modulus is constant and has the value 
which follows from 


20 


5G 3G 

— + 1 — — - 

Ef E 


(3.12) 


Et = 0.05725E when v = 0.3 (3-13) 

The bar is so dimensioned that the critical compressive stress for 
torsional buckling, in the absence of initial eccentricity (i.e. of 
initial twist), is by elasticity theory 


a = 2a* 


(3.1M 



16 


NACA TM 1592 


This is, at the same time, the value obtained by incremental theory. 
Deformation theory gives, however 


a = 1.05a* 


( 3 . 15 ) 


This result is calculated in supplement B. (in ref. 7 this value 
appears only in a plot.) Supplement B shows at the same time that 
E sec (ff) with 5 given by (5.15) is equal to 


E rp . = O. 56 OE when v = 0.3 


(3.16) 


With the inclusion of an initial eccentricity in the calculation, 
the definition of -what is considered as buckling load is naturally some- 
what arbitrary. In the case in point the maximum load which the bar is 
able to sustain (even without eccentricity the critical load is at the 
same time the maximum load) can be taken for this purpose. 

As Initial angle of twist per length 2b of the section 0.000111, 
0.00111, and 0.0111 rad. Is successively used. The divers results are 
found in table 3*9 (computed with the data of figure 9 from ref. J) . 

With this it is proved that, when proceeding from incremental theory, 
even a slight Initial eccentricity is sufficient to produce a considerable 
displacement of the computed critical load and that the value of it is 
comparatively little dependent on the Initial eccentricity (naturally 
only for an eccentricity already different from zero). The thus com- 
puted critical loads are close to those obtained by deformation theory 
with and without initial twist, for there is very little difference 
between them. 

The unusual feature of the Onat and Drucker load case is that In 
spite of the presence of Initial twist, the state of stress remains 
uniform during the application of the load. In other problems on plates 
‘which have an initial eccentricity, this is not the case as a rule, and 
problems such as these are very little amenable to calculation. However, 
it Is to be expected that then specific changes in the critical buckling 
load (which then indicates the load at which, more or less suddenly, the 
deflection assumes a significant form) will occur. 

Whether these changes will then show the same degree of independence 
from the actual amount of initial eccentricity, and whether the data then 
are again in close agreement with deformation theory (without counting 
initial eccentricity) has, of course, not been established as yet. It 
is not known even in the particular load case of Onat and Drucker for 
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values of the tangent modulus different from those which they used, 
and which is rather at the lower limit (Et ~ 0.05E) or for other dimen- 
sions of the structure. 

In this bu ckl ing case the stress-strain relations themselves exhibit 
the marked effect of initial eccentricity. (Cf. table 3*6 at Ej. = 0.05E.) 

(In (3.13) Efc = 0.05725E.) 

At T = O.10Q 


E(d7/dT) F = 4.260 


hence substantially higher than that at t = 0 (2.600). 

At a r value only, slightly above t = O.IcTq the value obtained 
for E(d 7 /dr) F is equal to E(d 7 /dT) D at r = 0 and E sec (a 0 ) = 0-5E 
(in (3.16), E sec (a 0 ) = O. 56 OE). 

However, it already has been proved that for other values of Ej. 
and E sec (a 0 ) the quotient E(d 7 /dr) p increases slower at increasing 
r and reaches the value of E(d 7 /dT)jj at t = 0 not until later. 

Nevertheless, the buckling case of Onat and Drucker gives a good quali- 
tative picture of the influence of initial eccentricities. 

Bijlaard and Wiseman likewise point out (in ref. 46) that the approxi- 
mate coincidence of the results by incremental theory, when initial eccen- 
tricity is allowed for, with those by deformation theory (initial eccen- 
tricity disregarded) in the Onat and Drucker buckling case is not applicable 
to other relations. 

That, all the same, deformation theory is always in agreement with 
experiment, as Bijlaard and Wiseman claim (cf . section 4.1), must have 
been known to them from the fact that incremental theory is not correct, 
but that deformation theory is valid, at least at load changes such as 
occur in buckling cases and where the data of the theory of Batdorf and 
Budiansky are equal to that of a deformation theory (not that by Nadai, 
compare section 3.3). (See section 10.1 of ref. 3, for example.) 

Bijlaard and Wiseman conclude their discussion of reference 7 with: 
"Nevertheless, it is a welcome contribution for indicating that small 
divergences in the direction of incremental theory are cancelled by 
initial eccentricities." 
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3.6. COMPARISON OP DATA FOR. CRITICAL LOADS BY DEFORMATION 
AND INCREMENTAL THEORY 


Such, a comparison is shown in figure 3.4, taken from reference 1. 

The reproduction of four load cases is the same as discussed in sec- 
tion 2. 3*2.3 • 

The dashed curves represent deformation theory, the solid curves 
incremental theory, both allowing for the Shanley effect . The plot 
also shows to which load cases the respective curves refer. 

Figures 3-5 and 3*6, taken from reference 9j give the same data as 
figure 2.1. The variation of the buckling factor k is plotted verti- 
cally, the ratio a/b of the sides of the rectangular plate hinged 
along all edges, horizontally. The plate is loaded on the side of 
length b in compression in one direction. 

• Figure 3 -3 represents the curves by incremental theory and, according 
to Bleich, both computed without the Shanley effect. The curve by elas- 
ticity theory is also shown. All curves were computed with consideration 
to the exact elastic value for the transverse contraction coefficient _ 
v = 0.32. 


Because Hahdelman and Prager wanted ter compare their results with 
that of Ilyushin, they computed the data with v = 1/2 ._ The particular 
comparison Is shown in figure 3-6. 

Figures 3-7, taken from reference 30, and 3-8, taken In part from 
reference 30 and in part from reference 31 > represent the data of several 
theories for an infinitely long plate hinged along the edges compressed 
in length direction. The shape of figure 3*7 is that according to fig- 
ure 2.4, that of figure 3*8 according to figure 2.3* - - 

The results of the several methods based on deformation theory are 
not far from each other and near the curve according to Gerard's method. 
The inclusion or exclusion of the Shanley effect does not matter much. 
(Compare Ilyushin's curves in figs. 3*7 and 3-8 with that by Stowell in 
the same plots.) 

Incremental theory gives consistently much higher data than defor- 
mation theory or Gerard's method. Whether the Shanley effect is Included 
or not, does not matter much. Compare the Handelman . and_ Prager curves in 
figure 3.8 with those of Bijlaard (in ref. 31) by incremental theory, 
computed with Shanley effect accounted for. 
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it-. EXPERIMENTAL DATA 


4.1. NACA TESTS 


Of the NACA tests for 1945-1946, -while the present report was 
written, references 32, 33 j 34,.35 j and 36 were available. These experi- 
ments deserve special attention since Bijlaard in reference 23 sees in 
the results of these experiments a confirmation of deformation theory. 

They are compression tests on H, Z, and channel sections. The 
section forms used, as described in reference 32, are found In figure 4.1 
those of references 33 to 36 in figure 4.2. The dimensions of the col- 
umns were so chosen that column buckling did not occur, but only local 
instability. In references 33 to 36 it is stated that: "The lengths of 

the columns were selected so as to obtain whenever possible a desirable 
three-half wave buckling pattern." In the description of the tests 
nothing is ever mentioned about the buckling pattern in flanges or webs, 
only the critical load is indicated. The cross-sectional distortion of 
the sections according to figure 4.1 is given once. (See fig. 4.3.) 

The selfsame photograph of an H-section under compression published 
in references 31 to 36 shows buckled flanges. The buckling pattern Is 
so laid out that It consists of three-half waves. The photograph shows 
no sign of web buckling. About the condition of the critical load the 
following statement is given: "In the local instability test, measure- 

ments were taken of the cross-sectional distortion, and the critical 
stress was determined as the stress at the point near the top of the 
knee of the stress -distort ion curve where a marked increase in distor- 
tion first occurred with small increase of stress." It is well known 
that such methods of determining the buckling load contain a large meas- 
ure of arbitrariness. (See ref. 48, for example.) 

In the tests of references 32 to 36 stress-strain curves are deter- 
mined on test pieces taken on different areas of the cross section. In 
each of the cited publications, either an average diagram is given which 
indicates the amount of scattering over all test pieces, or 5 to 10 dia- 
grams, which again were the averaged result of a number of test specimens 
from the same section bar. What seems to be about the extreme limit of 
the thus published stress-strain curves is represented In figure 4.4 to 
figure 4.7. It is seen then that Gerard's method (section 3.2) always 
gives somewhat higher results. The writers themselves apply no theory. 
Figures 4.4 to 4.7 simply contain curves which indicate the average loca- 
tion Of the test points. 

Reference 38 again discusses the results of the NACA tests of refer- 
ence 33, as far as the H-sections are concerned. Figure 4.8 represents 
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experimental points in a diagram of the form of figure 2.3. [The calcu- 
lation of the critical load "by elasticity theory is carried out according 
to reference 40. Theoretical curves are also shown. The theoretical 
critical loads in the plastic range are likewise computed according to 
the principles of reference 40 and with the aid of Stowell's report 
(ref. 21). 

A summary of references 32 to 3 6 is given in reference 39* From 
this reference .39 it is apparent that references 44 and 45 belong in 
the test series of references 32 to 36 . The last-named reports were not 
available, however, to the writer at the time this report was written. 

Bijlaard (in r&f. 22) gives the diagram on figure 4.9, in. which the 
curves for two load cases computed by deformation theory are shown. ^ 

Then, he quotes "NACA tests 1945-194 6." Of what tests precisely, is not 
mentioned. It is certain that the material was "avianal," but which of 
the materials used by Pride and Heimerl was meant by it? The quotes are 
not immediately traceable to the work of Pride and Heimerl. Neither is 
the stress-strain curve for "avional," which was copied by Bijlaard with- 
out the scatter of Heimerl and Pride. (A line is drawn through the eight 
points of figure 4.9 in reference 22, but that is no theoretical curve; 
it merely indicates an averaged position of the intended points.) Appar- 
ently Bijlaard did not compute the critical load by elasticity theory 
again, but simply followed the report of Pride and Heimerl. 

If, In the NACA tests, the flanges really have started to buckle, the 
values of the critical load should lie between load cases I and II, at 
least on the premise that the buckling flanges in the tests are still 
considered as infinitely long (three half-wave waves occur, according 
to fig. 4.9)* But even Gerard's simple rule itself yields results which 
are fairly accurate, for the experimental points are located near to the 
stress-strain curve of figure 4.9* 

Pride and Heimerl (ref. 30) also made somd conpression tests on a 
long, square, hollow, seamless box. The individual plates could then 
be considered as simply supported at the edges. 

Figure 4.10 represents test points for seven groups of test pieces, 
each group referring to test pieces from one sectional bar. Besides the 
test points, the curveB of Gerard, i.e. stress-strain curves, are plotted. 
There still is some scattering between the stress-strain curves of the 
different groups. And even within a group of test pieces from one bar 
the scatter is fairly wide. (Cf. the limits indicated on curve C.) 

4The support conditions of the ends infinitely far from each other is 
of no significance. For case II, in which °° waves occur, this Is Instantly 
clear. In case I a half-wave occurs f6r hinged ends; for clamped ends, 
a whole wave. But in this case the critical compressive load is indepen- 
dent of the wave length, if the wave length Is great. 
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The same measurements are repeated in figure 4.11, but after cor- 
rection (according to ref. 47), hence refer to material with a constant 
stress-strain curve. The theoretical curves are those of figure 3*5 • 
The agreement with Stowell’ s theory is good. Gerard’s secant-modulus 
method should also give acceptable results. 


4.2. GERARD’S EXPERIMENTS 


Gerard (ref. 19) likewise made some tests on Z and channel sec- 
tions and established a satisfactory proof of his secant modulus method. 
The chart in which the real critical load is plotted against the calcu- 
lation by elasticity theory obtained by secant modulus method, passes 
nicely along the test points. (See fig. 4.12 and fig. 4.13.) The cal- 
culation of the critical load by elasticity theory was made according 
to reference 4l. One noticeable feature of figure 4.13 is that the 
agreement between theory and experiment in the elastic range is not 
good,i which raises some doubt as to whether all experiments are not at 
the upper limit. The true critical load is determined from the differ- 
ence in strain at the opposite sides of a flange. At first the differ- 
ence seems to be practically zero, but at incipient buckling this dif- 
ference grows quickly at slightly increasing mean compressive stress 
(fig. 4.14). Figures 4.12 and 4.13 also show the curves according to 
Lundquist (ref. IS) and Laiighaar (ref. 17) (tangent modulus method). 

Other tests by Gerard are reported in reference 42 and discussed 
also by Stowell (in ref. 43). The tests refer to a shearing teBt on an 
infinitely long rectangular plate. The long edges of the plate are 
rigidly clamped. Reference 43 gives a diagram in which the critical 
shearing stress is plotted against the critical shearing strain computed 
by elasticity method, according to Stowell’ s theory and Gerard’s secant 
modulus method (fig. 4.15). Gerard’s curve forms the shearing stress- 
shearing strain curve of the material. Stowell' s curve can also be 
applied to the case of hinged long sides. (Naturally, that of Gerard 
is too.) 

This case was also treated by Stowell, but there was only a very 
sma.ll difference from the case of clamped edges, which cannot be taken 
into account in the calculation. 

The agreement between the experimental evidence and the curve by 
Gerard's method iB not as good in figure 4.15 as in figures 4.12 
and 4 . 13.5 

5 " • , - - . . 

Translator's note: Figure 4.15, taken from reference 43, appears 

to be in error and is misleading, Gerard's method actually gives lower 
stresses than Stowell' s method. (See MCA TN 3184.) 
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4.3* Kollbrunner 1 s Experiments 

In reference 8 Kollbrunner investigated rectangular plates of dif- 
ferent length/width ratios for several, load cases (fig. 4.l6). The 
edges of length b were assumed to be hinge supported, those of 
length a ifree, hinge supported or clamped, as illustrated in fig- 
ure 4.16. jThe tests were made on perfectly simple plates; thus, for 
example, the plate for case II (fig. 4.l4) was not the flange of an 
equilateral angle section loaded in compression nor the plate of case IV 
a side of a square box beam, such as the other investigators used. 

Figure 4.17 shows a portion of kollbrunner 1 s setup. The lateral, 
edges of the plates are enclosed in small steel tabs, which safeguarded 
the hinged support of the edges as illustrated in figure 4.18(a). The 
set screws in the steel tabs insured clamping (fig. 4.18(b)). 

Kollbrunner checked his experimental results against his theoretical 
solutions obtained by the differential equation (3*7) which is a kind of 
mean of Bleich's and Chwalla's relations. 

The theoretical result for the critical compressive stress 0 is 



with 

a plate length 

b plate width 

h plate thickness 

a by (3*4) and (3*5) 

m number of half-waves of buckling pattern 


(4.1) 
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The constants p and q. in the different load cases (fig. 4.l6) 
have the following values: 

Case I p = 0 q. = 0 


II 

p = 0.425 

<1 = 0 

III 

p = 0.570 

QL = 0.125 

IV 

(M 

II 

Pi 

H 

II 

o* 

V 

P = 2.5 

q. = 5 

VI 

p = 2.270 

q. = 2.450 


Coefficient K in (4.1) is a correction factor defined by experi- 
ment. The neccessity for the use of this coefficient is attributable 
to the fact that the assumption of hinged plate support on the short 
side where they are compressed is not perfect, since small bending 
moments can still be taken up here. 


The correction factor K for case I is not given. For case II a 
factor K = 1.5 should be applicable in the elastic as well as in the 
plastic range. For cases III to VI a factor K = 1.2 is applicable, 
but in the plastic range only. 


By (4.1) it can be deduced that the critical compressive stress 5 
for the infinitely long plate is 


a 



1 + Va 
2 



(4.2) 


as given by Kollbrunner himself. 

Coefficient K does not appear any longer in (4.2), and actually, 
the manner of support of ends loaded infinitely far from each other are 
of no interest in load cases II to VI. Compare case II with case I 
mentioned in the footnote of section 4.1. 


Kollbrunner does not report why the factor K = 1.5 for case I 
holds In the elastic and the plastic range, whereas for cases III t.o VI 
the factor K, which, moreover, has another value, is valid in the plastic 
range only. All the same the use of the correction factor K seems a 
rather arbitrary interference in the experimental check of a theory which 
itself has as yet no solid basis. (Cf. section 3*2.) On top of that, it 
should be remembered that the significance of K decreases considerably 
when long plates are involved. The figures 4.19, 4.20, and 4.21, of the 
same pattern as figure 2.1, represent some test data with the "theoretical" 
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curves, i.e. curves obtained with correction factor K, together with 
the curves of Bleich and Chwalla and Kollbrunner, respectively, of refer- 
ences 13 and 15. These graphs still show points which are up to 10 per- 
cent lower than KoUbrunner 1 s "theoretical" curves, which Kollbrunner, 
indeed, confirms in reference 11. 

To define the instant of buckling, a straight bar is reflected in 
the plate surface . According to Kollbrunner, "the instant of buckling 
was momentarily plainly visible by a typical jump of the reflected bar 
from its straight shape into the wave pattern." 

After the appearance of reference 20, Kollbrunner becomes a sup- 
porter of Ilyushin's theory, to which he subscribes, as stated before, 
with great enthusiasm in reference 11. Ilyushin's theory gives up to 
10 percent higher values than Kollbrunner' s theory, and the latter 
believes that if the correct (elastic value) is used instead of v = l/2 
in Ilyushin's theory, both theories would be in good agreement. 

Figures 4.22 and 4.23* of the same form as figure 2.2, represent 
some experimental results, based on the theoretical curves of KoUbrunner 
and Ilyushin. 


3. DISCUSSION 


The deformation theory of plastic strain applied to the calculation 
of the critical buckling load of plates gives results which are in better 
agreement with experiment than incremental theory. 

Onat and Drucker give an explanation for the failure of incremental 
theory in computing critical buckling loads for a special case by assuming 
an initial eccentricity. 

Some writers consider deformation theory as "perfectly correct" and 
are willing to see more in it than a practical calculating rule. It 
never has been proved even for general cases that incremental theory 
applied to a plate with initial eccentricities gives exactly the same 
results as deformation theory applied to the plate, without counting 
with the initial eccentricity. 

The statement of some writers that their or someone else's experi- 
ments are in "excellent agreement"' with their own theory needs to be 
critically examined. 

At such a pronouncement the question often arises: 

1. How much scatter was there in the stress-strain curves? (Cf . the 
scattering in figs. 4.4 to 4.7, for example.) 
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2. Exactly how is the buckling limit defined? Compare, for example, 
references 32 to 36 , where "the critical load is considered as reached, 
when the distortion of the cross section begins to show a rapid increase 
with slightly increasing stress." 

3. Li ttle scatter in the measured critical state of stress is no 
indication of whether there are disturbing secondary effects in the 
experiment or initial eccentricity. Groups of identical plates, loaded 
fairly far in the plastic range, almost all exhibit the same critical 
buckling load, but fairly far in the plastic range the tangent modu- 
lus Efc is us ually small, and the stresses which can occur even before 
buckling can then still not vary very much. However, should not the 
scatter in distortions on reaching the critical state be much greater? 

The inf o rma tion about these problems is generally very vague. Con- 
sidering reference 30 as example (already discussed in section ^-. 1 ), 
the following may be stated: 

To 1. The scattering in the stress-strain curves for one of the 
square boxes used as the test pieces seems to be considerable (fig. ^-. 10 , 
curve C.) Between the mean curves for the different boxes there are 
differences of the same kind. 

To 2. As examples of the determination of the buckling load two 
suitable diagrams are given in reference 30 (fig- 5-1)- But they refer 
to cases in which the buckling load was, at the same time, approximately 
the maxi mum load. Diagrams for cases in which the maximum load was con- 
siderably above the buckling load were not given. In this connection 
the question arises as to how it was possible that on many tubes the 
buckling stress and the mean stress at maximum load were almost identi- 
cal, but far below the 0.2 strain limit, as seen from table I of refer- 
ence 30 (it- 2.8 ksi, ^ 3.2 ksi, and 61.4- ksi). 

Lacking further information, this is a very remarkable result. The 
inescapable conclusion was that for a better insight into the problem 
we had to carry out experiments of our own. 


6. CONCLUSIONS 


1. It may be assumed that the plastic distortion increases even 
during buckling. Unloading according to elasticity theory is not 
required (Shanley effect). 

2. The incremental theory of plastic deformation is really the one 
which should be applied to the calculation of critical buckling loads 
of plates. 
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3. Initial eccentricities are bound to have a marked displacement 
effect on the result. As a rule, however, it is practically imp ossible 
to allow for such eccentricities in the derivation. 

Deformation theory applied without regard to initial eccentrici- 
ties gives practical results. 

5. Of all the writers starting from deformation theory, Bijlaard 
gives the best treatment. He accounts for the Shanley effect and intro- 
duces the correct value of the transverse contraction coefficient. 

6 . Uncertainties in the stress-strain curve at incipient buckli ng , 
likewise the scattering in the experiment, often show agreement between 
theory and experiment that seems better than it act ually is. 

7. The difference in theoretical results for the critical load on 
the basis of deformation theory with or without Shanley effect and with 
or without the correct value of v is not or, not much, greater than the 
scattering that is. to be observed as the result of many causes ringi ng 
the experiment. 

8 . Less explored is the critical distortion at which buckling occurs. 

• 9 . Gerard's method, which introduces the secant modulus, i.e. 

Young's modulus in the formulas for the buckling stress in the elastic 
range, always gives, it is true, a higher critical buckling stress than 
deformation theory, although it gives, nevertheless, a practical cal- 
culating rule. 

10. A better insight Into the problem Involved. may perhaps be obtain- 
able by personal experimental investigation. 


Translated by J. Vanier 
National Advisory Committee 
for Aeronautics 
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SUPPLEMENT A 


THE DIFFERENTIAL QUOTIENT d 7 /dr OF INITIALLY SIMPLY 
COMPRESSED MATERIAL AS FUNCTION OF THE LATER 
ON ADDED SHEARING STRESS t 

A.l. INTRODUCTION 

f 

The material is such that above the elastic l imit and at increasing 
load in a compression test the differential quotient do/de = E^ exhibits 

a constant value. The material is first considered as subjected to a 
compressive stress a = 0 q, after which at constant a , a shearing 

stress t is added. The secant modulus on reaching compressive 
stress (Tq is E 6ec (0Q), hence, also a constant in the subsequent 

calculations . 

A. 2. BY FLOW THEORY 

The formulas (ll.l6) and (11. 17) of reference 3 give the stress- 
strain relations for the case that only normal stresses 0 (positive 
if 0 is a compressive stress) and shearing stresses t are applied 

= -^ + | P 0 dd 2 (A.l) 

d 7 = 2 ^ 1 + dT + 2pr dip (A. 2) 

£> 

e the normal strain (positive when e is a compression). 

J 2 = | a 2 + t 2 (A. 3) 

and p is a function of Jg. First to be investigated is the function 
of Jg on the basis of constant tangent modulus E^ in the compres- 
sion test, where, because t = 0 
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and by (A.l) and (A. 3 ) 


or 


T _ 1 -2 

J 2 " 5 


de = £2 + £ pcr 2 da = de 

T x T TJ 


E 


®t 


(A.10 


_(A. 5 ) 


and with (A. 4 ) 




P = T^fe- - |) 


(A.6) 


4 J 2 Wt E> 

Formula (A.6), where 0 = 0 Q and (A.6) are substituted in (A. 2 ), gives 

i-2 

'o T 5t£ 

which is identical with (3«9) and at t = 0 changes to (3. 11). 


d/ _ 2(1 + v) 9 t 2 /I 1\ 

dT ” E 0 rt 2 + 5T 2 VFb " E / 


(A. 7 ) 


A. 3. BY DEFORMATION THEORY 


Formulas (H.l^) and ( 11 . 15 ) of reference 3 give relations between 
stresses 0 and ' t end strains s and 7 


e = SL + £ pa 
E 3 


(A.8) 


7 = - t Ji l t + Pt 

E 

with P a function of 3 0 with J 0 according to (A. 3 ). 


(A. 9 ) 


Now assume that the equation of the stress-strain relation for the 
compressive test ‘in the plastic range is 
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This, substituted in (A. 8), gives 



and since (A. 4) is valid again for this test , P becomes 


(A. 10) 



(A. 11) 


Formula (A.3) in 'which 


cr = ( and (A. 11) are substituted in (A. 9) 


gives 


T - + V> T + 

' E 



3s-t 

\ + 5 t2 


(A. 3-2) 


The constant a can be eliminated by ascertaining that after 
applying the compressive stress Oq according to (A.IO) 



is valid, which in (A. 12 .) giveB 


2(1 + v) 


E 


t + 3 


(st ' S 


T + 


3ao T 


ft 


+ 3t £ 


sec 


( * 0 ) 


*t) 


- (A.15) 


The differential quotient beccmes 


d7 

2(1 + 

v) A 3 

. 3 + . 

3gq 

ft _ gr 2 V 1 . l\" 

dT 

E 

Et 

E 

h 2 + 5t2 

a Q 2 + 3 T ^ / /\ E sec^ a o) ^ 


(A. 14) 
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hence is identical with (3*8) and at t = 0 becomes (3»10) and when 
E sec( a 0) = chaiiges to (3. 11). 
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SUPPLEMENT ~R 

DETERMINATION BY DEFORMATION THEORY OF THE CRITICAL BUCKLING 
LOAD OF THE ONAT AND DRUCKER LOAD CASE IN THE ABSENCE 
OF INITIAL ECCENTRICITY AND PREDICTION OF THE 
SECANT MODULUS FOR THIS LOAD 


In the Onat and Drucker buckling case the yield point is a*. The 
critical compressive load by elasticity theory is 5 e = 20*. Above the 

yield point the tangential modulus E^ is assumed constant, and has the 
value t'hat follows from 


20 =^ +1 -¥ (b.U 

hence Et = 0.05E for v = 1/2 and Ej. = O.O57E for v = 0.3. 

Above the yield point the relation between compressive stress and 
compression is 


a 0 * a* 
Et E 


(B.2) 


and the secant modulus is 


E 


sec 


a 

Et 


a* a* 

Et E 


(B.3) 


At small increase of shearing stress t the differential quotient 
by deformation theory is identical with (3.10) 


sec 


= d 7 = 2v - 1 
dT E 


3/_ct_ _ O* + £*\ 
cr\pt ^ E J 


(B.4) 


For the special load case in question the critical value of a by 
incremental theory has the value according to elasticity theory, hence 
0 = 5 e = 2c*, as long as no initial deflections exist. 
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The critical value for a by deformation theory becomes 


G, 


2G C 


cr = rsec 5 = ^§ec a * 


(B.3) 


In view of (A.l 8 ) and (A.19), the value of 5 is the solution for 
a from the equation 


£_ _ G( 2v - 1 + 3_ 
a 2\ E Et 


_ 3 _ sl + 1 Si 
Et a Ed 


a = 


2 + 2. 


3 

I 


2 v - 1 . 3 


1 ■ r ^ 


+ 4 - ^ + 1 


-¥ 


+ 1 


(B. 6 ) 


and with (B.l) 


a = 1 . 05 a 


When the compressive stress reaches this value , the secant modulus 
is, by (B.3) and (B.l) 


E sec = °- 525 E for v = l /2 


E 


sec 


= O. 56 OE for v = 0.3 
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TABLE 3.1. 

COMPARISON OF THE DIFFERENTIAL QUOTIENT dy/dr BY INCREMENTAL 
THEORY (F) AND DEFORMATION THEORY (D) OF MATERIAL UNDER 
INITIAL CCMPRESSIVE STRESS a Q IN THE PIASTIC RANGE 


Ej. = 0.8E 


T 

*<&), 

• 

4 

ll) 

E sec( a 0) = E 

E sec( CT 0) - °-9E 

0 

2.600 

2.600 

2.933 

. OOIOq 

2.600 

2.600 

2.933 

.003a 0 

2.600 

2.600 

2.933 

• OIOq 

2.600 

2.600 

2.933 

.03^0 

2.602 

2.603 

2.935 

0 

b 

H 

• 

2.622 

2.632 

2.951 

• 3^o 

2.759 

2.826 

3.059 













58 


NACA TM 1392 


TABLE 3.2. 


CCMPARISON OF THE DIFFERENTIAL QUOTIENT d 7 /dT BY INCREMENTAL 
THEORY (F) AND DEFORMATION THEORY (D) OF MATERIAL UNDER 
INITIAL COMPRESS IVE STRESS a Q IN THE PLASTIC RANGE 


Ej. = 0 . 6E 


T 

! 

■©l 


i 7 N 

1t/ d 

® MO (”0) = B 

E sec( 0 o) = 0.9E 

0 

2.600 

2.600 ■ 

2.955 

• OOIcTq 

2.600 

2.600 

2.955 

. 0030q 

2.600 

2.600 

2.955 

• Olcro 

2.601 

2.601 

2.954 

•03oq 

2.605 

2.608 

2.940 

• 

H 

Q 

O 

2.658 

2.687 

5.005 

O 

t> 

• 

5.025 

5.203 

5.435 
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TABLE 3-3. 

COMPARISON OF THE DIFFERENTIAL QUOTIENT d.y/dr BY INCREMENTAL 
THEORY (F) AND DEFORMATION THEORY (D) OF MATERIAL UNDER 
INITIAL COMPRESSIVE STRESS a Q IN THE PLASTIC RANGE 


% = 0.4E 


T 

■ 

E (^2S 
\ dT /D 

E sec( CT 0) = E 

EgecC^O) = 0.9E 

E sec( CT 0) = 0.5E 

0 

2.600 

2.600 

2.933 

5.600 

• 0010*0 

2.600 

2.600 

2.933 

5.600 

. 003ctq 

2.600 

2.600 

2.933 

5.600 

«01oq 

2.601 

2.602 

2.935 

5.601 

• O30q 

2.612 

2.618 

2.950 

5.606 

.10Q 

2.731 

2.795 

3.114 

5.665 

•3cr 0 

3.557 

3.956 

4.189 

6.052 
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TABLE 3.4. 

COMPARISON OF THE DIFFERENTIAL QUOTIENT (17/dT BY INCREMENTAL 
THEORY (F) AND DEFORMATION THEORY (D) OF MATERIAL UNDER 
INITIAL COMPRESSIVE STRESS a Q IN THE PLASTIC RANGE 


E t = 0.2E 


T 

E (t4 


E ee=(°0) = S 

S 8ec ( a o) = 0.9E 

E S ec(*0) “ 0.5E 

0 

2.600 

2.600 

2.955 

5.600 

• OOIcJq 

2.600 

2.600 

2.955 

5.600 

• 003cr 0 

2.600 

2.600 

2.955 

5.600 

• OIcTq 

2.6o4 

2.605 

2.958 

5.6o4 

• 03cTq 

2.632 

2.648 

2.980 

5.636 

1 * la o 

2.949 

i 5.120 

5.439 

5-990 

•3°0 

5.151 

6.216 

6.448 

8.312 
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TABLE 3.5. 

COMPARISON OB’ TBE DIFFERENTIAL QUOTIENT CLy/ir BY INCREMENTAL THEORY (F) AND 
DEFORMATION THEORY (D) OF MATERIAL UNDER INITIAL COMPRESSIVE 
STRESS a 0 IN THE PLASTIC RANGE 


Et = O.IE 

T 

E&l) 

E 

(SA 

E Bec( a 0) = E 

E sec(°o) ** O.SE 

EsecC^o) - 0 . 5 E 

®sec(^o) - 0 . 15 E 

0 

2.600 

2.600 

2.933 

5.600 

19.600 

. OOlffg 

2.600 

2.600 

2.933 

5.600 

19.600 

.003c q 

2.602 

2.600 

2.933 

5.600 

19.600 

.oia 0 

2.608 

2.611 

2.944 

5.610 

19.604 

.° 3 c 0 

2.673 

2.708 

3.040 

5.696 

19.640 

• lcr 0 

3.386 

3-770 

4.089 

6.640 

20.033 

• 3a 0 

8. 340 

10.735 

IO.968 

12.831 

22.613 
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TABLE 3 . 6 . 

COMPARISON GF TEE DIFFERENTIAL QUOTIENT dy/dr BY INCREMENTAL THEORY (F) AND 
DEFORMATION THEORY (D) OF MATERIAL UNDER INITIAL COMPRESSIVE 
STRESS ff 0 IN THE PLASTIC RANGE 





Ej. = 0 . 05 E 

- - 

T 

e (S) f 

e(J21) 


E se c( a 0) =E 

E aec( tT 0) =- 0 . 9 ® 

Bsec^O) = 0 . 5 E 

E sec(°o) = 0 . 15 E 


0 

2.600 

2.600 

2.935 

5.600 

19.601 

• 001a 0 

2.600 

2.600 

2.933 

5.600 

19.601 

• 003cfo 

2.601 

2.600 

2.933 

5.600 

19.601 

• OIcTq 

2.617 

2.623 

2.956 

5.622 

19.617 

• 05 c iq 

2.753 

2.828 

3.160 

5.816 

19.761 

.IOq 

It-. 260 

5.070 

5.389 

7 - 94-0 

21.334 

« 

V* 

Q 

O 

14.718 

19 . 77 ^ 

20.007 

21.870 

31.653 
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TABLE 3.7. 

COMPARISON OF THE DIFFERENTIAL QUOTIENT dy/dT BY INCREMENTAL THEORY (F) AND 
DEFORMATION THEORY (D) OF MATERIAL UNDER INITIAL COMPRESSIVE 
STRESS ff Q IN THE PLASTIC RANGE 


Efc =• 0.02E 


T 

e (S) f 


E sec( a o) - E 

E sec( ff 0 ) = 0 . 9 E 

E sec( cr o) “ 0.5E 

E aec( ff 0 ) - 0 . 15 E 

0 

2.600 

2.600 

2.933 

5.600 

19.601 

• OOIcXq 

2.600 

2.600 

2.933 

5.600 

19.601 

• 003^q 

2.601 

2.600 

2.933 

5.600 

19.601 

• olcr o 

2.6H 

2.659 

2.992 

5-658 

19.653 

• 03<Jq 

2.996 

3.188 

3.520 

6.176 

20.121 

* la o 

6.882 

8.970 

9.289 

11.810 

25.231 

•3cr 0 

33.891 

*1-6.891 

*(-7.121- 

18. 987 

58.770 
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TAB IS 3.8. 

COMPARISON OF THE DIFFERENTIAL QUOTIENT d^/dr BY INCREMENTAL THECRY (F) AND 
DEFORMATION THEORY (D) CF MATERIAL UNDER INITIAL COMPRESSIVE 
STRESS c 0 IN THE PLASTIC RANGE 


E,. = O.OIE 

. T 

E (S) f 

E 

fez) . 

Vdr/ D 

E sec( a O) = B 

W^O) = 0 . 9 E 

E S ecfe 0 ) = 0 . 5 E 

Esecfeo) » 0 . 15 E 

0 

2.600 

2.600 

2.933 

5.600 

19.601 

.OOIcTq 

2.601 

2.600 

2.933 

5.600 

19.601 

• OOJCTo 

2.608 

2.600 

2.933 

5.600 

19.601 

•°l ff O 

2.689 

2.719 

3.052 

5.718 

19.713 

■ 03^0 

3.4oo 

3.788 

4.120 

6.776 

20.721 

0 

t> 

H 

• 

11.251 

15.470 

15.789 

18.340 

31.734 

• 3 ®o 

65.74-2 

92.086 

92.319 

94.182 

103.96 
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TABLE 3.9. 


RESULTS OF 0NA1 AND TRUCKER BUCKLING CASE 
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Figure 2.1.- Critical compressive stress a x of rectangular plate 
loaded under compression in one direction, all edges hinge 
supported; length of sides a and b; width-depth ratio (b/h) is 
constant. The sides of length b are the loaded sides, 
m = half-waves. 
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Figure 2.2- Method of plotting used by Kollbrunner. Critical com- 
pressive load a x of an infinitely long plate with hinged edges 
under compression in the direction of the length. 
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Q008 
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f 


Figure 2*3.- NACA method No. 1. Vertical: critical stress cr x 

of an infinitely long plate of 14S-T6 aluminum with hinged edges 
under compression in length direction. Horizontal: critical 
strain F x>e by elasticity theory. 




Figure 2.4.- NACA method No. 2. Load case as in figure 2.2. Verti 
cal: _ factor tj by which a x e should be multiplied according to 
elasticity .theory to yield actual critical stress a x . Horizontal: 
actual stress a 71 
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Figure 3.4- Buckling factor k plotted against 7 X (or shear 7 ) for 

four different cases of 24S-T aluminum, according to Besse ling 
(ref. 1). 
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k=<T x 


I2(t-y 2 )b 2 

7T^Eh2 



0 


b 


Figure 3.5.- Buckling factor k for simply supported rectangular 
plate v =0.32 hinge supported on all edges under compression. 
The sides of length b are loaded, width/thickness ratio according 
to b 2 /h 2 = 1000. m = half-waves at buckling (taken from ref. 9). 
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figure 3.6.- Buckling factor k for simply supported rectang ular 
plate v = 0.50, all edges hinge supported, loaded in compression. 
The sides of length b are loaded; width/thickness ratio according 

to b2/b2 = 1000. m = half-waves at buckling (taken from ref. 9). 
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Critical stress <T % , ksi 


Figure 3.7.- Factor ti plotted against critical stress a x for an 

infinitely long plate of 14S-T6 aluminum plate, with hinge -supported 
edges under compression in length direction (taken from ref. 30). 
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Figure 3.8.- Load case in figure 3.4. Critical stress a x plotted 
against critical strain ? x>e computed by elasticity theory (taken 
from refs. 30 and 31). 
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Figure 4.2.- Sections used in tests described in references 33 to 36. 



Figure 4.3.- Cross-sectional distortion during buckling in tests 
described in reference 32. 
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Scatter of stress-strain 



Critical stress calculated according to. 
elasticity theory, ksi 


Figure 4.4. - Experimental results of compression tests on sections 
described in reference 33; material, 75S-T aluminum. 



Critical compressive stress, ksi 


60 


NACA TM 1392 



Critical compressive stress 
by elasticity theory, ksi 


Figure 4.5.- Experimental results of compression tests on sections 
described in reference 34; material, 24S-T aluminum. 
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Scatter of stress- strain 



Figure 4.6.- Experimental results of compression tests on sections 
described in reference 35; material, R 303-T extruded al uminum. 
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0 20 40 60 80 100 120 140 160 

Critical compressive stress 
by elasticity theory, ksi 


Figure 4.7.- Experimental results of compression tests on sections of 
14S-T extruded aluminum (from ref. 36). 



Critical compression according 
to elasticity theory 

Figure 4.8*- Experimental results of compression tests of reference 
31, as concerns H-sections; compared with theoretical results of 
Stowell (taken from ref. 38; for bp, b^, tp, and tw> see fig. 4.2). 
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I Infinitely long, rectangular plate 
comprssed in direction of length*, 
one unloaded edge hinge sup- 
ported, one unloaded edge free. 

31 Infinitely long, rectangular plate 
compressed in the length 
direction; one unloaded edge 
clamped, one unloaded edge free. 
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r or fastening 
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ends see 
footnote, 
section 4.1. 
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Critical compression according 
to elasticity theory 


Figure 4.9.- Two buckling cases calculated by Bijlaard, compared with 
NACA test on sections (taken from ref. 22). 



Critical compression according to elasticity theory 

£ x ,e = 

I2(l-v 2 )b 2 

Figure 4.10,- Critical stress a ^ ol infinitely long 1 14S-T aluminum 
alloy plate, tinge -supported edges, compressed in length direction. 

The seven curves refer to seven groups of test pieces, each group 
being manufactured from the same box -beam section. The curves 
are those of Gerard, hence stress-strain curves (from ref. 30). 
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Critical stress according to elasticity theory, ksi 

Figure 4.12.- Results of compression tests on 24S-T aluminum sections 

by Gerard (taken from ref. 19). 



Critical stress, ksi 
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Figure 4.13.- Results of compression tests on sections of 75S-T 
aluminum alloy by Gerard (taken from ref. 19). 
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Critical shearing stress according 
to elasticity theory 


Figure 4.15.- Critical shearing stress t_ of an infinitely long plate 
loaded in shear; theoretical curves for hinged and clamped edge 
support. Tests with long sides clamped. 
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Loaded edges are regarded 
as hinge- supported 
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Figure 4.16.- Kollbrunner load cases on avional m alloys. 
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(a) Hinge type 


(b) Fixed- clamping type 



Figure 4.18. - Steel sides used by Kollbrunner. 
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A = Ilyushin 
B = Kollbrunner 


Figure 4.23.- Critical buckling stress for load case V of figure 4.16 
on infinitely long plate according to Kollbrunner and Ilyushin. 
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V 




Figure 5.1.- Determination of critical load and maximum load in the 
compression tests of reference 30. (Results of these tests given in 
figs. 4.10 and 4.11.) 


NACA - Langley Field, Va. 


